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Abstract. We consider transformations of deterministic and random signals 
governed by simple dynamical mappings. It is shown that the resulting signal 
can be a random process described in terms of fractal distributions and fractal 
domain integrals. In typical cases a steady state satisfies a dilatation equation, 
relating an unknown function f{x) to f{Kx) (for example, f{x) = g{x) f (kx)) . 
We discuss simple linear models as well as nonlinear systems with chaotic be- 
havior including dissipative circuits with delayed feedback. 
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One of the most interesting and significant questions in nonlinear dynamics is 
the problem of noise influence on chaotic motion. Traditional deterministic models 
provide a fair description of averaged motion only for stable trajectory systems. 
In case of unstable motion random fluctuations may grow, so that non-stochastic 
descriptions fail which asks for deeper understanding and investigation of noisy 
dynamics. 

In this work we consider an evolution of simple linear and nonlinear noisy discrete 
dynamical systems. We concentrate on the important peculiarities of the stochastic 
evolution: self- similarity symmetry and fractal distribution generation. 

Consider a stochastic mapping: 



in real d-dimensional space. Assuming being a sequence of uncorrelated random 
variables, we can write the Kolmogorov-Chapman equation, which for large TV — > oo 
reduces to the equation for the stationary (asymptotic) distribution Pst that can 
be written in the form: 



The corresponding equation for the distribution function Fourier transform 



1. Introduction 



(1) 



Xn+1 = Cat + F{Xn) 



(2) 



dYdZ Pfi{X - Y) 5{Y - F{Z)) Pst{Z). 




is given by 

(3) *.t(t/) = *fl(f/) J dVaiU,V) -^stiV), 

(4) aiU.V) = (2^)-'^y dYeM^{Y,V)~t{FiY),U)}, 
where the scalar product takes the form: {X, U) — xiUi + • • • + XdUd- 

2. Easy cases: linear deterministic and noisy mappings 

We can consider a mapping ([T]) as a transformation of a d-dimensional signal in a 
circuit with dissipation. It means that the mapping is contractive and its Jacobian 
satisfies the condition | det{dFi/ dxj}\ < 1. A steady state solution corresponds to 
the stationary regime of deterministic or stochastic motion. As a preliminary we 
briefly discuss several simple cases. 

• Linear deterministic dissipative mapping. Let us suppose that a signal X 
subject to a dissipative contraction: X kX , k < 1, and mixed with an external 
signal: X ^ A + X . In this case we have 

F{X)^A + fiX, 0<K<1, 
Pfi{X)^S{X), 

and Eqs.([2]), ([3]) take the form of dilatation (scaling) equations: 

P,t{X)^K-^Pst {k-\X~A)) , 
^st{U) = exp{-i{A, U)} *st(Kf/). 

Utilizing a normalization ^'st(O) = 1 and using an iteration procedure we find the 
solution: 

^st{U) = f[ exp{-tK^{A, U)} = exp{~t{A, U) k^} = exp i-t^A^X , 

7=0 7=0 ■' 

Pst{X)^6{X~A/{l~K)), 

which describes coherent accumulation in a circuit. 

• Linear stochastic dissipative mapping with a Gaussian noise term. 

Assuming 

F{X) = A + kX, 
Pa{X)=P^{X,R) = (nR)-"/^ exp{-{X,X)/R} , 
and solving the equations 

Pst{X) - K-'' J dVPgiX - Y, R)Pst {k-\Y - A)) , 

^st{U) = exp{-ii?(C/, U) - i{A, U)}^stiKU) 
using the same procedure we obtain the solution 

*st(C/) = n '^^P {-^'^^(A U)~-^Rk'-'{U,U)^^ exp - i i? 



PstiX) =Ps{X- A/ {I - «), R/{1 - ^)) 
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A„ Ai K A3 X* A4 As 

(a) (b) 

Figure 1. (a) The scheme of branching for X = 6. (b) The tree 
for the thhd generation of a pre-fractaL 

It can be viewed as a resuh of two independent additive processes - for both co- 
herent and random components of the signal. 

• Linear stochastic dissipative mapping with Gaussian and Kubo- Andersen 
stochastic terms. Let us consider a model with a mixed stochastic term: = 

+ ^/v"^- I'^i t^is ca.se: 

F{X) ^ fiX, 

L-l 

Pfl = Pg * Pka, Pka{X) = J2 PkS{X - Xk), 
where * denotes the convolution. Now Eqs.(l2]),(l3]) take the form: 

(5) PstiX) = - / dYP,{X -Y,R)Y, PkPsi 

k^o 

(6) 4'st(C/) = exp{--i?([/,C/)}{^PfeexpH(Xfc,C/))}«'st(Kf/). 

The solution of Eq. ([6]) by the iteration procedure leads to a distribution with fractal 
properties: 

(7) = exp{-ii?|^} \{[Y^pueM~^^^ {Xk.U))] 

7=0 k=0 

It is convenient to represent this solution through integral over the fractal domain, 
or multifractal integral (MFI). The MFI concept is discussed in the next section. 

3. Multifractal integrals 

Consider the fractals with L-branch multiplication (the case of dichotomous 
branching was considered in 1 ) . Having assigned the values 

A, e(0,l), OeeAo<Ai<---Al_i = 1 

and the contraction ratio k < 1, we choose the law of branching, shown in Figure 
1, and define the multifractal integral of f{x), with x E [0, 1], as 

(8) / f{x)d^,ix\K, 8) = lim a„, an = i"" V eW/(AW), 



(existence and equality conditions for the limits in ([5]) are formulated below). The 
summation in ^ is made over all possible signature vectors s„ = [siP, ■ ■ ■ ,si"-*) 
with components Sn' € {Ag, Ai, . . . , Ax,_i} , where A^, = A^, — A^. The arguments 
of the function / are defined as follows 

(9) AW = A, + (l-«)(s„,h„(K)), 

where h„(K) = (1 ^) and (a, b) — aibi-\ |-a„6„, and the quantities 

[si 

Qn are functions of s„. Choosing normalization condition f7„ = 1 for / = 1 we 
find that 

(10) Veri = LeW, 



fc=0 



for = s„_i © (Afe) (where the © denotes a "concatenation" of vectors). In this 
section we assume that the conditions 



L-1 



(11) ew = n^(^»')' E^(^,) = L, 



k=0 



(12) Im^/;(Afe) = 0, Rei^(Afc)>0, 

hold. In this case one can interpret pk — ip{Ak)/L as probabilities; the correspond- 
ing multifractal measures induced by the multiplicative Besikovitch process were 
considered in Setting Xi — {1 — k)(A* + Aj), we represent the solution of Eq.([5]) 
(the Fourier transform of (O) through the MFI 



Pst{x) = J p^ix-x)d^i{x\K,e). 



c 

The following existence conditions hold for MFI: 

Theorem 1. If a function f{x), x G [0,1], is continuous and conditions (jlip . 
P2p hold, the sequences of the quantities (t„ in ([5]) converge to a limit equal to MFI. 

Proof. For certain n, the summation in ([5]) goes over the points of the n-th 
generation pre-fractal, which is the set of L" points with coordinates A^f' G (0, 1). 
By arranging the pre-fractals on parallel lines and connecting the points according 
to the rule s„ ^n+i — Snffi(Afc), fc = 0, . . . L— 1, we get a graph that represents the 
formation of the Cantor set for n — > oo. Let the graph contains pre-fractals of first 
n + k generations. One can see that every point of the n-th generation pre-fractal 
is a parent vertex for a cluster containing a tree of the subsequent generations. The 
corresponding rule of branching is given by 

ALt:' = AW + All = (1 - .)." (r,, h(.)) ; 

-k"A, < < k"(1- A,). 

As fc — > oo, every "small" cluster converges to a Cantor fractal, which is similar to 
the full fractal. 

Using definition ([8]) and conditions (fTTjl . we obtain 

(13) \(7n+k - <yn+k' \ < If n+fc — Cr n\ + Wn+k' ~ CTn], 

(14) a„+, - a„ = Ei-"0L^1 E^^'Qi-' (/(^«' + ^nl) - /(Ak^l)) . 
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Taking into account that \Qn^\ < 1, we get an estimate 

(15) - a„| < 5]L-"eW 5]L-'=eM [/(A^ + A^,) - /(A^) 

s r 

Let f{x) be continuous on the closed interval [0, 1] (thus, uniformly continuous): 

Vxe[0,l] Ve>0 3(5e (Va : |a| < (5,) : \f{x + a)-f{x)\<e/2. 
Using the condition 

AW - ^"A., AW + «"(! - A.)] c [aW - S,, aW + 

we can select TVe for every (5e. Employing the relations pTjl and P^ - P^ . we obtain: 

Ve > 3N^ yn > : |cr„+fc - (Jn+k' \ < e, 

showing that {(t„} is a Cauchy sequence and thus converges. I 

We find an alternative representation for the MFI using the singular distributions 
approach. Introducing singular functions 

(16) 9Jl„(x|K,e) = L-"y eW5(a;-AW), SW(a;|K, 6) = lim Srn„(a;|K, 9), 

s 

we have: 

1 1 
f{x)dfi{x\K,e) ^ I f{x)Tl{x\K,e)dx = lim [ f{x)Tln{x\K,e)dx. 



The Fourier transform of the limit in (|16p is a solution of a dilatation equation 
OT(u;|k, 9) = m{Lu\K, 9) M(kw|k, 9), 



L-l 



fn(u;|K, 9) = L-^ J2 {-it^(l - k)(A, + Ak)} 

and is equal to 

oo 



L-l 



7=0 L fe=0 

Accordingly, the dilatation equation for the limit function (jl6p reads: 



oc 

an(a:|/t, 9) = J m{x ~ y\K, 9) m {K~^y\K, 9) dy, 



L-l 



OT(x|k, 9) = ^ *(Afe) (5 (x - (1 - k)(A, + Afe)) 



It is possible to define the MFI as a limit of a sequence of definite integrals. Such 
definition is more illustrative but in fact is equivalent to ((8|)- (fT2)) . This approach 
was developed in [T] for L — 2. Here we present this formalism briefly. 

Define the MFI as 

(17) J f{x)df,{x\K, 9) = Jim a*, a* = 2- ^ 9^ <f (aW)^,„ , 



where 

-ifix^s = 6-' / f{t)dt 

Jx-5/2 

and aL"'^ is defined by ^ with i = 2, A* = i, Ao = 0, Ai = 1. For k < ^, a* is 

computed by integration over a system of nonintersecting segments £„ = IJ^. C 

[0, 1], where In^ = [Xn^ — ^k", aIT' + ^k"]- The set £„ is obtained from the unit 
segment by the standard Cantor procedure repeated n times: one first removes from 
the unit segment a middle part of length 1 — 2k, then the same fraction is removed 
from the middle of each of the remaining segments, and so forth. The measure is 

[si . fsl . . 

constructed by assigning to each segment In a weight factor 0„ appearing in the 
summation in (|17p . 

The helpful intuitive picture of the pre-fractal domain £„ and the cluster set 
C ~ linin^oo -Cn embedded in [0,1] is no longer valid for ^ < k < 1, since the 

[si 

segments In intersect in that case. This difficulty may be resolved by replacing 
the segments with rectangles and "spreading" them along the second coordinate. 
Define two-dimensional MFI of f{x, y), x e [0, 1], y G [0, 1] as follows 

(18) J I fix,y)df,ix,y\K,,Ky,e) = Ji^m 2-"^eW (/ (aW, AW))^^ . 
We assume that Aq]i is determined by the expression that follows from ([9]) by 

[si 

replacing k with Kq, a = x,y, and that the average is taken over the region dn = 
Ixl (8) 4" J lai = [ AqIi - ^K^, AqIi + ]. The rectangular cells dn^ do not overlap if 
at least one of the contraction coefficients Ka, ct = x,y, is not larger than i. This 
condition is satisfied, in particular, for k^, = k, < k < 1, and Ky — ^, and then 

J J fix)d^i{x,y\K,^,e) ^ J f{x)d^i{x\K,&). 
V c 

If < ^ for a = X (for a — y), the projection of I? = lim„^oo Us onto the 
X (the y) axis is a standard Cantor set that has the similarity dimension Da = 
— (log2Ka)~^. As it was shown in [5], the Hausdorff-Besicovitch dimension D of 
the whole set V is equal toZ? = 2/(l — log2 k). Note that for k = ^ we have D — 1; 
in this case 2? is a diagonal of the unit square, and (flSl) reduces to a definite integral 
on the diagonal. 

4. The nonlinear Ikeda mapping 

Consider a circuit with a nonlinear element and a delayed feedback [H [31 21 [51 [SI 
[3 [8]. Assuming (,^{t) being a random process with zero mean value {{^^{t))) = 0, 
we can write a stochastic difference equation 

(19) X{t) =. F{X{t - Tdei)) + e^(i - Tdei), 

where Tdci is the delay time. For X{t) being a two-dimensional (complex) signal, 
the equation (jl9p can be rewritten in the form of mapping (U]), where Xj^ = X{ta + 
NTdei), d = C(io + (A^-l- l)Tdei), < to < Tdei (in this case {X, U) - Re {XU*)). 
Suppose that the external noise is the Ornstein-Uhlenbeck random process with 
small correlation time Tcor Tdei (the "rapid" Gaussian noise). Assuming the 
transformation of a slowly varying complex valued amplitude X{t) as a phase shift 
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0— > if) + X\X\^ + 00 and the dissipative contraction \X\ k\X\, k < 1, we find the 
nonhnear term in ([T]) , (|19p in the form 

(20) F(X) = l + KXexp(iA|X|Vi6io) 

(the Ikeda mapping [3l|4]). As shown in [5l[6l[7], the random phase approximation 
can be used in the intence phase mixing limit A ^ 1, which allows to replace the 
delta-function in ^ by 

5(2) (y _ p(^z)) ™> ({S^^\Y - 1 - AtZe'*)))^, 

where {{■■■)) 4, denotes an average over 0. Then 

PstiX) = j dY P^{X - Y,R/{1~ n^)) P,^{Y). 

The radially symmetric distribution Pch describing noise caused by dynamic chao- 
tization [SJ [HI [7] can be written in the form 



(21) Peh(^) - (27r«:2)-i / I3dl3 S(/3) Jq {Pk^^ \X - 1|) 

Jo 

CXD 

(22) S(/3) = Umi^P). S(/3) = n UPn''). 

k=0 

The integrand in (|2T|) contains a factor S(/3) which is a solution of the dilata- 
tion equation ([22]) . In order to simplify (I2ip . ((22|) . we approximate the Bessel 
function by a function with separated "oscillatory" and monotonic parts Jq{x) sa 
x(a;) [v^ cos [x — jt^)] ■ In this approximation 

/eh(Q)«^|e-*-/4y' <|.(^VQ+?^,q) dl,{x\K,Q) + C.C. 

- 00 

'i>(/3,Q) = ^f3x{f3VQ) n ' = 2"/2exp [-iz^(s„, h„(l))] , 

fc=0 

where $ is the Fourier transform w.r.t. first argument of the monotonic function 
Notice that quantities 0„ fail to satisfy the conditions p2|l . Thus the Theorem 
1 does not hold and we have to look for an alternative condition of convergence. 

5. Extended class of multifractal integrals 

Let us consider the definition ([8]) of the MFI without the condition (fT2l) . Now 
we assume that ^'(Afc) are complex- valued. 

Theorem 2. Let the function f{x), x € [0, 1], be dijferentiable arbitrary many 
times, also let 3M Wx Wl € N : |/'-'-*(2:)| < M, and let the conditions (|lip hold. 
Then the sequences an in ^ converge to a limit equal to MFI. 

Proof. Taking into account the formula ^ and introducing the power expansion 

00 

/(AW) = E ^/('HA,)(1 - k)' (s„, h„(.^))' , 
1=0 

we obtain 

00 

(23) k„| < ^ ^|/(')(A,)| (1 - kY eW (s„, h„(«))' 



The use of the generating function for the powers of the scalar product exp{^(s„, h„(K))} 
and the Leibnitz formula allows us to perform the transformations: 

L-"5:eW (s„, h„(.))' = |,(L-"Eel^' e^^-^^"^'^))) 



L-l 



«=0 



k=0 



«=0 



L-1 



(24) =E n('^™^'-'"'^[^"'E*(^fc)^r^ 

m j=l fc=0 

The left summation in goes over all non-negative integers that satisfy the 
condition X]"=o ™i ~ ^'^'^ ("^i ■ ■ ■ "'•n)' = (™i + ' ' ' + '7i„)!/(mi! • • • m„!) are the 
multinomial coefficient. The inequality |Afc| = |Afc — A*| < 1 implies |A™| < 1, 
m € N. This allows us to obtain the upper estimates 

L-l L-l 



fc=0 



< max 

l/3fcl<l 



A;=0 



Gi < G, 



where G = max{Gi, 1}. Notice that the number of multipliers L J2k=o ^(Afc)^™ 



with m 7^ in (I24|) does not exceed the least of the numbers n and This allows 
us to write a sequence of upper estimates: 

(25) 

n L-l 

L-"E©"' h„(^))'| < E (mi, ■ . .m„)! J] ('*™^^'"'^ E *(Afc)Ar- 



fc=0 



(26) 



<G' E {rni,...mny. < G'(1 - k)" 

As result, we see from ([23]), ([25]) and ^ that 

oo ^ oo ^ 

(27) |a„| < E ^l/^"(A*)| G' < Af E ^yG' = Me^. 

Taking into account the relations ([121)- (US]), ([23 -([211), wc find: 

- |E ^""0"' E (/(^n' + ^n.k) - 

s r 

CX) ^ 

^ |E^(E^""0i^^'/^'^(^i^^)) {E^-'Qt'i^'Ay) 

1=1 ' s 

°° 1 < I 

^E7f E^""0L^'/''^(^i^^) '^"'(i--)' E^"'4^'(r^'h.(.))' 



oo ^ 



We see that 

Ve > 37V, Vn > TV, : \an+k - (Jn+k' \ < 2Me^ (e^"" - 1) < e, 
which implies convergence of the Cauchy sequence {an}- H 



6. Random processes transformations in a circuit with dissipation 

and delayed feedback. 



Consider a circuit described by the equations p^ . (|20p in assumption that 
T cor ^ Tdei [SJ [H [71 [5] . This situation is more difficuh for analysis as one should 
take into account correlations between the variables — ^(tg + {N + l)Tdci) with 
different N. Now we are forced to seek the stationary solution of the generalized 
Kolmogorov- Chapman equation for the multi-time distribution functions: 

(28) P,t((X.),e«+i) = J clYd4P,tiiYs),Co) 

n 

In (|28p we use a short-hand notation 

Pjv((^s),e) = P((Xo[7Vrdel], Xi[7VTdel + Ti], • • • ,X„[iVTdel + r„]), elA^Tdel]), 

where Xi, i — ...,n, are the signal amplitudes at NT del + Ti, respectively, and 
^ denotes the amplitude of at NT ^d- We also require = tq < n < • • • < 
T„ < Tn+i = Tdei and use a notation = (IAq dAi dA2 .... The noise term is 
assumed to be the Markovian random process with arbitrary r cor and the transition 
density function w(^,^', [t]) (the Ornstein-Uhlenbeck and Kubo- Andersen random 
process were treated in [HIIZIIH])- The equation for the distribution function Fourier 
transform reads 

(29) Vst((C/.),r!„+i) = J dVdili;,t{{Vs),no) 

n 

X Y\_<^iUp,Vp) H{np+i,np - Up,[Tp+i -Tp]), 

p=0 

where H{fl, fl' , [r]) is the Fourier transform of w{£,, [r]) and the function a{U, V) 
is defined by (|4]). This function may be written in the form 

(T = a + Aa, a {U, V) = {2tt \V\y' e-'^-'^J {\V\ - k \U\) , 

where the first term a corresponds to the random phase approximation and de- 
scribes the case of chaotic motion with intense phase mixing [71 [S]. This means 
that in the limit A ^ cx) in Ij20p one can replace a by a (respective approximate 
solution of ([28]) . ([29]) are labeled by superscript (0)). The exact and approximate 
forms of ([^5)) can be presented in the operator form: 



random phase approx. (o) (o) 



In case of the Ornstein-Uhlenbeck random process [8] the explicit formula for K 

takes the form: 

(31) 

K/((C/,), C/„+i) = [ n e-'«'^^''] $(([/.), C/„+i) (( /((Kf/«e'^0, E 

p=0 \\ 9=0 

„ n+1 n+1 a 2 0^ 

mu.), [/„,,) = exp{-| E IE J^M (l - ^) }. 

p^l q—p 'J'p-i 
9 




It is evident that the operator K produces a scaling transformation (dilatation), 
because the right-hand side of (j3ip contains the functions with scaled arguments 
kUs- Thus the right equation in (|30p is a generalized dilatation equation. The 
equation (|28p may be viewed as an equation with the "smaU" operator eS. This 
enables to obtain the solution in the form of a series expansion. Taking into account 
the normalization conditions -0 st ((0) , 0) — V'^'st ((0) i 0) = 1? finds 

oo oo 

vi? m) , n) = jhn^ K™ 1, Vst = 0it°^ + E { (E s} Vi? • 

This approach was developed in J8] where the criterion of convergence was exam- 
ined. 
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